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ABSTRACT 

Motivated by recent detections by the XMM and Chandra satellites of X-ray line emission from hot, 
luminous stars, we present synthetic line profiles for X-rays emitted within parameterized models of a 
hot-star wind. The X-ray line emission is taken to occur at a sharply defined co- moving- frame resonance 
wavelength, which is Doppler-shifted by a stellar wind outflow parameterized by a 'beta' velocity law, 
v(r) — i>oo(l — R*/r)P . Above some initial onset radius R Q for X-ray emission, the radial variation of the 



emission filling factor is assumed to decline as a power-law in radius, f(r) 



The computed emission 



profiles also account for continuum absorption within the wind, with the overall strength characterized by 
a cumulative optical depth r* . In terms of a wavelength shift from line-center scaled in units of the wind 
terminal speed , we present normalized X-ray line profiles for various combinations of the parameters 
/3, r*, q and R Q , and including also the effect of instrumental and/or macroturbulent broadening as 
characterized by a Gaussian with a parameterized width a. We discuss the implications for interpreting 
observed hot-star X-ray spectra, with emphasis on signatures for discriminating between "coronal" and 
"wind-shock" scenarios. In particular, we note that in profiles observed so far the substantial amount of 
emission longward of line center will be difficult to reconcile with the expected attenuation by the wind 
and stellar core in either a wind-shock or coronal model. 

Subject headings: radiative transfer — shock waves — stars: early-type — stars:mass-loss — X-rays: 
stars — X-rays: line-profiles 



1. INTRODUCTION 

Over the past two decades, orbiting X-ray observatories 
from Einstein to ASCA have shown that both hot and cool 
stars are moderately strong sources of soft X-rays at ener- 
gies betweenm 0.1 and 2 keV. For cooler, late- type stars, 
the scaling of X-ray luminosity with, e.g., stellar rotation, 
suggests a solar-type origin, with surface magnetic loops 
confining hot, X-ray-emitting material within nearly static 
circumstellar coronae. For hotter, early-type stars, the X- 
ray luminosity has no discernible dependence on rotation, 
but instead scales nearly linearly with bolometric luminos- 
ity, L x « 10 _7 Lbo1 (Long & White 1980; Pallavicini et al. 
1981; Chlebowski, Harnden, & Sciortino 1989), suggest- 
ing a different emission mechanism, perhaps related to the 
high-density, radiatively driven stellar wind observed from 
such stars. 

Nonetheless, an ongoing question has remained whether 
such hot-star X-ray emission originates from coronal-like 
processes occuring in nearly static regions near the stellar 
surface (Cassinelli & Olson 1979; Waldron 1984; Cohen, 
Cassinelli, & MacFarlane 1997), or rather in the highly su- 
personic wind outflow, perhaps from shocks arising from 
intrinsic instabilities in the radiative driving mechanism 
(Lucy & White 1980; Lucy 1982; Owocki, Castor, & Ry- 
bicki 1988; Feldmeier et al. 1997; Feldmeier, Puis, & Paul- 
drach 1997; Owocki & Cohen 1999). Based on the ob- 
served strong winds in hot stars (e.g., Snow & Morton 
1976), a key discriminant has been the apparent general 
lack of soft X-ray absorption, which, along with optical 



and UV spectral constraints, implies limits to the extent, 
temperature, and fractional contribution of the total X- 
ray output from a base corona (Cassinelli, Olson, & Stalio 
1978; Nordsieck, Cassinelli, & Anderson 1981; Cassinelli & 
Swank 1983; Baade & Lucy 1987; MacFarlane et al. 1993; 
Cohen et al. 1996). This limited inferred absorption is 
thus generally seen as evidence for emission arising from 
lower-density regions at large radii, in general agreement 
with a wind-shock model. 

In the context of extreme ultraviolet spectroscopy, Mac- 
Farlane et al. (1991) presented models of emission line pro- 
files, from both a coronal source and from an expanding 
shell, affected by transfer through a stellar wind. These 
authors demonstrated that emission line profiles are a 
key discriminant of the source location. However, due to 
the small throughput and modest resolution of the Ex- 
treme Ultraviolet Explorer (EUVE) spectrometers, com- 
bined with the extreme opacity of the interstellar medium 
in this bandpass, emission lines from only one hot star (e 
CMa) were observed with EUVE. The only high signal-to- 
noise line seen in this EUV spectrum was the He II Lya 
line, which was observed to be broadened by several hun- 
dred km s _1 , lending credence to models that postulate a 
wind-origin for the hot plasma in OB stars (Cassinelli et 
al. 1995). 

It is only recently, however, that orbiting X-ray observa- 
tories like XMM and Chandra have been able to provide 
the first detections of spectrally resolved X-ray emission 
lines from hot, luminous stars (Kahn et al. 2001; Schulz 
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et al. 2001; Waldron & Cassinelli 2001; Cassinelli et al. 
2001). Through the broadened profiles of these resolved 
X-ray emission lines, these data provide a new, key diag- 
nostic for inferring whether the X-ray emission originates 
from a nearly static corona, or from an expanding stellar 
wind outflow. A key goal of the present paper is to pro- 
vide a firm basis for interpreting such X-ray emission line 
spectra. 

Specifically, we describe below the method (§2) and re- 
sults (§3) for computations of synthetic line profiles aris- 
ing from parameterized forms of X-ray emission in an ex- 
panding stellar wind. In comparison to recent work by 
Ignace (2001), who derived analytic forms for X-ray emis- 
sion profiles for the case of constant expansion velocity, 
our analysis takes into account a variation of velocity v(r) 
in radius r, as parameterized by the so-called 'beta-law' 
form v(r) = Voo(l — R*/r)@ , where R* is the stellar radius 
and Voq is the wind terminal speed. As in Ignace (2001), 
the overall level of attenuation is characterized by an inte- 
gral optical depth r* , and the X-ray emission filling factor 
is assumed to decline outward as a power law in radius, 
f(r) ~ r~ q . However, we also allow for this emission to 
begin at a specified onset radius R Q that can be set at 
or above the surface radius R*. We first present scaled 
emission profiles for various combinations of the four pa- 
rameters /?, t*, q, and R /R* (§3.1). For selected cases 
intended to represent roughly the competing models of a 
"base coronal" emission vs. instability-generated "wind- 
shock" emission, we next compute synthetic profiles that 
are broadened by the estimated instrumental resolution for 
emission lines observed by XMM and Chandra (§3.2). Fi- 
nally, we discuss (§4) how various features of the synthetic 
profiles are influenced by properties of the emitting mate- 
rial, and outline prospects for future work in interpreting 
hot-star X-ray spectra. 

2. COMPUTING X-RAY LINE EMISSION FROM A 
SPHERICAL STELLAR WIND 

Extended monitoring observations by X-ray missions 
like ROSAT (Berghofer & Schmitt 1994) have shown that 
for most early-type stars the X-ray emission is remarkably 
constant, with little or no evidence of rotational modu- 
lation, and with overall variations at or below the 1% 
level. This suggests that the hot, X-ray emitting material 
must be well-distributed within the circumstellar volume, 
whether in a narrow, near-surface corona, or as numer- 
ous (n > (1/0. 01) 2 = 10 4 ), small elements of shock- heated 
gas embedded in the much cooler, UV- line- driven, ambi- 
ent stellar wind. Here we model any such complex struc- 
ture in terms of a spherically symmetric, smooth distri- 
bution of X-ray line emission, with emissivity rj\(r, /i) at 
an observer's wavelength A along direction cosine p from 
a radius r. The resulting X-ray luminosity spectrum L\ is 
given by integrals of the emission over direction and radius, 
attenuated by continuum absorption within the medium, 
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2.1. Optical Depth Integration 

The absorption optical depth t[/i, r] can be most easily 
evaluated by converting to ray coordinates z = fir and 



p = T\J\ — /i 2 , and then integrating over distance z for 
each ray with impact parameter p, 



t\p, z] 



Kp[r f ] dz f , 



(2) 



where p[r') is the mass density at radius r' = \j p 2 + z' 2 . 
We assume here that the X-ray absorption cross-section- 
per-unit-mass, k, is a fixed constant that is independent of 
wavelength A within the line and radius r within the wind. 
For emission at a given r along direction cosine /i, the op- 
tical depth attenuation in eqn. (1) can then be obtained 
through the simple substitution, 

r[p,r]=t ^l- p 2 r, pr . (3) 

For simplicity, we consider here a steady-state wind with 
velocity given by the 'beta-law' form w(r) = v(r)/voo = 
(1 — R*/r)P . Taking account also of the occult ation by the 
stellar core of radius i?*, we then write 

; p<R* & z< ^ R* 2 - p 2 

R*dz' n . , A . 

(4) 
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Here we have used the mass loss rate, M = Aitpvr 2 , 
to define a characteristic wind optical depth, r* = 
kM/AttVoqR*. Note that in a wind with (3 = and thus a 
constant velocity v = v^, the radial (p = 0) optical depth 
at radius r would be given simply by t[0, r] = r*R*/r. 
Thus in such a contant- velocity wind, r* would be the ra- 
dial optical depth at the surface radius i?*, while R\ = 
t*R* would be the radius of unit radial optical depth. 

For non-integer /?, the optical depth integral (4) must 
generally be done by numerical quadrature. But for gen- 
eral integer values, e.g. (3 = 0, 1, 2, or 3, it is still possible 
to obtain closed-form analytic expressions for this optical 
depth, though the expressions can be quite cumbersome 
to write down explicitly (see also Ignace (2001)). In the 
present work, these analytic forms are derived and eval- 
uated using the mathematical analysis software, Mathe- 
matica. To provide a specific example, we cite here the 
analytic result for the canonical case of (3 = 1, for which 
the integral in eqn. (4) becomes 



arctan 



z' jr' 



arctan 



(5) 

This expression actually applies under the full general con- 
ditions cited in eqn. (4), though some care must be taken 
in its evaluation for the cases with p < 1. 

2.2. Line- Emission Integration 

Let us next consider the line emissivity, which we take 
here to be of the form 

Vx (r, M ) = Cp 2 (r)f(r) 5 [A - A (l - fxv{r)/c)} , (6) 

where C is an emission normalization constant, A G is the 
rest wavelength for the line-emission, and f(r) represents 
a volume filling-factor for X-ray emission. Note that the 
radial dependence of this filling factor may reflect varia- 
tions in the temperature distribution of the plasma, as well 
as variations in the total emission measure. 
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The form (6) simply assumes the emission is an arbitrar- 
ily sharp function centered on the local, Doppler-shifted, 
comoving-frame wavelength A G (1— fiv(r)/c), and so ignores 
all broadening (e.g., natural, Stark, thermal, turbulent). 
For Doppler shifts associated with either thermal or small- 
scale "micro-turbulent" motions, the delta-function should 
be replaced with a finite-width emission profile. But gen- 
erally, such microscopic motions can be expected to be no 
greater than the gas sound speed, which for X-ray emitting 
plasma of temperature T x ~ 10 7 K is of order a few times 
100 km s _1 . As such, the associated broadening should 
be relatively small compared to that expected from the 
directed motion of the wind outflow, which has a much 
larger characteristic speed > 1000 km s _1 . Moreover, 
as discussed in §3.2, for larger-scale, "macro-turbulent" 
motions, the effects can be taken into account a postori 
by incorporating them with the general smearing of the 
observed line profiles expected from the limited spectral 
resolution of the X-ray detector. 

For convenience, let us define a scaled wavelength mea- 
sured from line center in units of the wind terminal speed 
^oo, i.e. x = (A/A — l)c/voo. In terms of this scaled wave- 
length, the integrals (1) for the luminosity spectrum then 
become 

/l roo 
d/i / drr 2 f(r)p 2 (r)e- T ^ r U(x+/iw(r)) , 

(7) 

where we have used the convention L x dx = L\d\. In- 
tegrating the Dirac J-function over direction cosine /i, we 
then find that the wavelength-dependent X-ray luminosity 
can be evaluated from the single radial integral, 

(8) 

where H is the Heaviside (a.k.a. unit-step) function, with 
H(s) = 1 if s > 0, and zero otherwise. 



Figure 1 plots the optical depth r[—x/w(r),r] versus 
scaled wavelength x at selected radii r in models with var- 
ious (3. For each radius, there is a critical, positive wave- 
length x where the optical depth increases steeply, repre- 
senting the strong attenuation of the red-wing emission by 
the dense inner wind and stellar core. Such core occupa- 
tion makes the integrand in eqn. (9) effectively zero for 
all u where (1 — u)P\Jl — u 2 < x. For blue- wing wave- 
lengths x < 0, the termination of the curves indicates the 
maximum blue-shift allowed by the velocity at the relevant 
radius. 

Ignace (2001) has recently derived completely analytic 
results for the X-ray profiles assuming a power-law emis- 
sion factor f(r) ~ r~ q and constant outflow velocity, (3 = 
(implying w{r) = 1). As noted above, for the somewhat 
more general case of nonzero, but integer values of (3, the 
optical depth integral (4) can still be done analytically; 
but in general the radial integral (8) must be evaluated by 
numerical quadrature. For this numerical integration, it is 
convenient (following Ignace (2001)) to define an inverse 
radius coordinate u = i?*/r, yielding 

where u x = 1 — l^ 1 /^. 

3. RESULTS FOR NORMALIZED LINE PROFILES FROM 
POWER LAW X-RAY EMISSION 

For any chosen model for the spatial variation of emis- 
sion /(r), eqn. (9) can be readily evaluated through a 
single numerical integration. As a specific illustration, 
we now present results for the example case of an emis- 
sion factor with the power-law form f(r) = f*u q for radii 
r = R*/u > R = R*/u Ql and zero otherwise. This cutoff 
implies a slight redefinition to the upper integral bound in 
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eqn.(9), 



min u Q , 1 



|V/3 



(10) 



To focus on the general form of X-ray line emission, let 
us define a normalized emission profile, 

«x = — \- v (11) 

max x [L x \ 

which simply sets the profile maximum to have a value 
of unity. This effectively removes the dependence on the 
several free parameters (e.g. /*, M, etc.) that af- 
fect the overall emission normalization. The form of the 
resulting X-ray line-profile l x versus scaled wavelength x 
thus depends on four free parameters, viz. /?, t*, g, and 
R /R*. The effects of wind attenuation and other model 
parameters on the overall X-ray emission level have been 
addressed in a broad-band context by Owocki & Cohen 
(1999), but in focussing here on line profile shapes, we de- 
fer to future work any analysis of the overall emission line 
strength. 

3.1. Synthetic Emission Line Profiles for Various Model 

Parameters 

Figure 2 presents synthetic emission profiles for vari- 
ous combinations of these parameters, all plotted vs. the 
scaled wavelength x over the range — 1 < x < 1. For 
comparison, figure 3 shows the flat- topped profiles for cor- 
responding cases with optically thin emission, r* = 0.01. 
Some notable general properties and trends are as follows: 

• The profiles seem most sensitive to r* and R Q1 with 
/3 and q having a comparatively modest influence. 

• The emission peaks are always blueward of line cen- 
ter (i.e., at x < 0), with increasing blueshift for 
increasing r* and/or increasing R Q . 

• The overall line-width likewise increases with in- 
creasing r* and/or R Q . The narrowest lines oc- 
cur when both r* and R Q are near unity. Large 
r* with low R tend to give blue-shifted "peaked" 
forms, while large R Q tend to be blue-tilted "flat- 
top" forms. 

• The blue-to-red slope near line center is always neg- 
ative. For low r* it varies from very steep for low 
R to nearly flat for large R Q . For large r* it has a 
more intermediate steepness for all R Q . 

3.2. Broadened Profiles for Coronal vs. Wind- Shock 
Scenarios 

For this simple parameterized model, the profiles pre- 
sented above represent the expected intrinsic line emission, 
with the only broadening associated with the directed, 
large-scale wind outflow. In practice, observed line pro- 
files can be further broadened, both by the limited spectral 
resolution of the detector, as well as by Doppler shifting 
from additional, random motions of the emitting gas. As 
already noted in §2.2, for small-scale thermal or "micro- 
turbulent" motions, the associated broadening can be ex- 
pected to be no larger than the gas sound speed, which for 
T x ~ 10 T K plasma is order a few 100 km s _1 . Since this 
is much less than the Vqq > 1000 km s _1 characteristic of 



the directed wind outflow, it should have comparitatively 
minor effect on line-profiles arising from the outer wind. 
For the inner wind or a more static "coronal" models, 
larger-scale turbulent broadening, which might in princi- 
ple have a larger characeristic speed, can still be accounted 
for a postori, together with instrumental broadening aris- 
ing from the limited spectral resolution of the detectors 
used in the XMM and Chandra satellites. For simplicity, 
we assume here that both macro-turbulent and instrumen- 
tal broadening can be characterized by convolving a simple 
Gaussian kernel function over the intrinsic line-profile. As 
such, a single parameterized width <j, defined here in units 
of the wind terminal speed , represents the combined ef- 
fects of both any macroturbulent motion and intrumental 
broadening, with a = V 'tf urb + <r? natr . 

Let us now consider examples of such broadened line- 
profiles. To keep the parameter range tractable, we iden- 
tify specific parameter subsets to represent general con- 
ditions in the two distinct scenarios that have often been 
invoked in modeling hot-star X-ray emission, namely the 
"coronal" (e.g., Waldron 1984) vs. "wind-shock" scenarios 
(e.g., Feldmeier, Puis, & Pauldrach 1997; Owocki & Co- 
hen 1999). The coronal scenario envisions X-ray emission 
as coming mostly from a relatively narrow, nearly static 
layer near the stellar surface; we represent this here with 
the parameter values q = 5, /? = 3, and R Q = R*. The 
wind-shock scenario generally envisions X-ray production 
arising from instability-generated shocks that form in the 
acceleration region r = 1 — 3R* of the star's line-driven 
wind, with eventual slow decay of the shocks in the outer 
wind; we represent this here with the parameter values 
R = 1.5i?*, (3 = 1, and q = 0.5. For both scenarios we 
again consider the effect of various optical depth parame- 
ters, T*. 

Figure 4 compares results for the coronal (left) vs. wind- 
shock (right) scenario. In characterizing results, we shall 
refer to profiles with small vs. large a as representing 
good/high vs. coarse/low resolution. Here some notable 
general properties and trends include: 

• At the best resolution (a = 0.1), the coronal mod- 
els are, as expected, notably narrower than the wind 
models. 

• The wind-shock models are generally more blue- 
shifted than the coronal models. 

• At the coarsest resolution, (a = 0.5), all profiles 
have similar width, but still generally differ in over- 
all blue-shift. 

• For the optically thin case, the emission is nearly 
symmetric about line center. For other cases, the 
emission peaks become more blue-shifted with in- 
creasing optical depth r*. (The coronal models 
with larger r* may not be relevant, however, since 
wind absorption would greatly reduce the overall 
strength of any observable emission.) 

The typical spectral resolution for Chandra is around 
A/AA ~ 1000 — 1500. For characteristic wind terminal 
speed of ~ 2000 — 3000, this implies a typical broaden- 
ing parameter of a ~ 0.1. For XMM, the spectral broad- 
ening is comparable, though slightly larger. Clearly, at 
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Fig. 2. — Scaled emission profile l x vs. the scaled wavelength x, plotted for various parameter combinations. In each box, profiles for 
optical depths t* = 1, 3, 5, and 10 are overplotted, with the various curves distinguished by having their peak values scaled respectively by 
factors of 1, 0.9, 0.8, and 0.7. The six rows represent, from the top to bottom, cases with initial emission radii of R /R* = 1, 1.1, 1.5, 2, 3, 
and 10. The central of the three columns represents a 'standard' case with (3 = 1 and q = 0, while the columns to each side illustrate cases 
which differ from this standard by setting (3 = 3 (left) or q = 3 (right). 
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Fig. 3. — Flat-topped emission profiles for the optically thin case r* = 0.01. The left, center, and right box correspond to the same wind 
parameters as in figure 2, but the five overplots in each box now represent cases with initial emission radii at R /R* = 1.1, 1.5, 2, 3, and 10, 
with the various curves now distinguished by having their peak values scaled respectively by factors of 1, 0.95, 0.9, 0.85, and 0.8. 



such high resolutions, it should be possible to place strong 
constraints on the site of X-ray emission and the degree of 
wind attenuation, as well as on the characteristic speeds 
of directed vs. turbulent motions of the emitting plasma. 

4. DISCUSSION 

The parameterized emission models here provide a good 
basis for interpreting X-ray line-profile observations made 
possible by the new generation of orbiting X-ray tele- 
scopes. While detailed, quantitative analyses must await 
the general availability of specific datasets, it is already 
clear from the few hot-star X-ray spectra published thus 
far that there are some noteworthy challenges for modeling 
the associated X-ray emission. 

For example, temperature- sensitive line-ratio diagnos- 
tics provide constraints on the magnitude of shocks in 
dynamical models. Generally instability generated wind- 
shock models tend to be dominated by moderate shock- 
strengths Av ~ 100 — 400 km s _1 that can only heat ma- 
terial up to a maximum of ~ 10 T K (Feldmeier et al. 1997). 
In those stars with much higher inferred temperatures, for 
example 6 1 Ori C with T > 5 x 10 7 K (Schulz et al. 2001), 
there is a need to produce stronger shocks. One possibilty 
is for wind outflow to be channeled by magnetic loops (Ba- 
bel & Montmerle 1997). In this case, wind material ac- 
celerated up opposite footpoints of a closed loop can be 
forced to collide near the loop top, in principle yielding 
shocks with velocity jumps of order Av ~ 1000 km s _1 , a 
substantial fraction of the wind terminal speed. 

Another key diagnostic stems from the forbidden-to- 
intercombination (f/i) ratios of helium-like emission lines, 



which can provide a strong constraint on the location 
of the X-ray emission. In the case of the solar corona, 
such ratios are nominally interpreted as density diagnos- 
tics (Gabriel & Jordan 1969). But because the strong UV 
radiation of hot-stars can generally dominate over elec- 
tron collisions in the destruction of forbidden line upper 
levels, in such stars the f/i provides instead a constraint on 
the strength of the UV intensity (Blumenthal, Drake, & 
Tucker 1972; Me we & Schrijver 1978). With appropriate 
atmospheric modeling of the relevant part of UV contin- 
uum, any inferred dilution of the radiation thus provides 
a diagnostic of the proximity of the X-ray emission to the 
stellar surface. For XMM observations of £ Pup (Kahn 
et al. 2001), initial interpretations suggest formation well 
away from the stellar surface (r > 2i?*), consistent with a 
wind-shock origin, with moreover a general trend toward 
larger radii for lower-energy lines, as would be expected 
from the larger optical depths and thus larger unit-optical- 
depth radii for such lines. For Chandra observations of 
£ Ori Waldron & Cassinelli (2001), the inferred formation 
radii are also generally away from the star (r ~ 2 — 10i?*), 
but with one stage (Si XIII) suggesting a formation very 
close to the surface (r < l.LR*). However, in this latter 
case, the relevant UV flux for destruction of the forbidden 
upper levels lies at energies above the Ly man-edge, and 
so any error in atmospheric modeling of the Lyman jump 
could still alter this inferred formation radius. Chandra 
observations of ( Pup (Cassinelli et al. 2001) show simi- 
lar results, with most inferred formation radii in the range 
r ^ 1.5 — 10i?*, with trends again consistent with inferred 
unit-optical-depth radii, but again with one stage (in this 
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Fig. 4. — As in figure 2, scaled emission profile l x vs. the scaled wavelength x, but now comparing results for a coronal (left) vs. wind-shock 
(right) scenario. The four overplots now represent profiles with broadening parameters set to a = 0, 0.1, 0.3 and 0.5, again with the overall 
emission scaled respectively by factor of 1, 0.9, 0.8, and 0.7. The five rows now represent (from the top) cases with optical depth parameters 
set respectively to r* = 0.01, 1, 3, 5, and 10. 
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case S XV) implying formation quite close to the surface 
(r < 1.25i?*). Of course, such fir constraints on source lo- 
cation must eventually be considered in conjunction with 
additional constraints, such as those provided by analysis 
of line profiles. 

Indeed, in the context of this paper, the most notable 
new diagnostic provided by recent X-ray observations lies 
in the form of the line profiles. For some initial results, 
these have some quite puzzling characteristics. For exam- 
ple, in the XMM spectrum of C Pup (Kahn et al. 2001), the 
lines all exhibit a similar, significant broadening, but with 
half-widths ~ 1000 — 1500 km s _1 that are still somewhat 
lower than the inferred wind terminal speed Voo > 2000 km 
s -1 . The typical profiles have a rounded, Gaussian form 
that is quite distinct from, e.g., the fiat-topped profiles 
expected for optically thin emission from a nearly con- 
stant velocity expansion. Moreover, the profiles show only 
modest asymmetry, with the peak intensity only slightly 
blue-shifted, by ~ — 400 km s _1 relative to line-center, 
and the red-wing only slightly flatter and more extended 
than the blue-wing. Recent somewhat higher spectral res- 
olution observations by Chandra do show the greater shift 
and asymmetry expected from optically thick wind emis- 
sion for C Pup (Cassinelli et al. 2001). But in Chandra ob- 
servations of the 09.7 lb star ( Ori, the profiles are again 
very broad (half- widths of ~ 1000 km s _1 ), but nearly 
symmetric (Waldron & Cassinelli 2001). 

Such profiles with substantial broadening, but little 
asymmetry and blueshift are difficult to reconcile with a 
simple wind-outflow model that includes attenuation by 
the wind and stellar core, both of which should preferen- 
tially attenuate X-rays from the back-side hemisphere, and 
thus in a wind or a corona with some outflow, preferen- 
tially reduce the red-side emission. In fact, from detailed 
modeling constrained by ROSAT observations, Hillier et 
al. (1993) have argued that that the wind of £ Pup should 
be optically thick to 1 keV X-rays at 10 i?*, with the ra- 
dius of optical depth unity being even greater for softer 
X-rays (20 R* at 0.6 keV; and still 2 R* at 2 keV). 

Even for a base corona with some overall outflow ve- 
locity and no overlying wind attenuation, occultation by 
the stellar core should still produce line profiles that are 
skewed toward the blue. A coronal model with little or no 
net outflow could explain the observed profile symmetry, 
but matching the observed line- widths would then require 
the coronal material to have extremely large "turbulent" 
motions. For example, line widths in the Chandra observa- 
tions of ( Ori (Waldron & Cassinelli 2001) require velocity 
dispersions of order 1000 km s _1 , and some of the lines ob- 
served by Chandra in 6 1 Ori C (Schulz et al. 2001) are even 
broader than this, yet apparently still nearly symmetric. 
If the lack of red-wing attenuation is interpreted to mean 
that red-shifted emission does not arise from wind out- 
flow from the backside hemisphere, then this suggest that 
material in front of the star must be flowing downward 
at velocities of 1000 km s _1 . Since this is greater than 
the surface escape speed, it cannot be the result of grav- 
itational infall, but instead would imply that some mech- 
anism, e.g., perhaps magnetic reconnection (Waldron & 
Cassinelli 2001), must actively accelerate roughly equally 
strong upflow and downflow with speeds of 1000 km s _1 ! 
Moreover, such random motions should inevitably lead to 



high-speed gas collisions and their associated shock heat- 
ing. For shocks with velocity jumps of order 1000 km s _1 
the associated post-shock temperature should be of order 
10 8 K, implying generally harder X-rays and higher ion- 
ization stages than are typically observed. In this context, 
these broad but symmetric profiles would indeed seem to 
represent a strong challenge for theoretical models. 

On the other hand, our synthesized line profiles from 
wind models do seem capable of producing relatively 
broad, yet symmetric and unshifted profiles, but only in 
the nearly complete absence of wind attenuation, and with 
emission originating far from the occultation of the stellar 
core. To test the viabilty of this alternative, observed line- 
profiles should be interpreted in conjunction with other 
diagnostics, like the f/i ratio, that can constrain the radii 
of X-ray emission, and with an accurate account of wind 
absorption, as computed through opacity calculations for 
the observationally inferred wind mass loss rate. 

Any wind shock model with significant wind attenuation 
is expected to be strongly blueshifted and quite asymmet- 
ric (see figure 4). Of all the line profiles discussed in the 
first group of XMM and Chandra papers, the observa- 
tions of Pup provide the best evidence for the blue-shifts 
and asymmetry expected for optically thick wind emission 
(Kahn et al. 2001; Cassinelli et al. 2001). For example, fig- 
ure 3 of Cassinelli et al. (2001) shows 5 out of 6 lines with 
blueward shifts and asymmetry similar to profiles calcu- 
lated here for wind-shock models with substantial optical 
depth, e.g. r* ~ 5. At least for this canonical O-type star, 
these profiles seems to provide strong evidence of a wind 
origin for most of the observed X-rays. 

In summary, systematic fitting of the line-profiles will 
be necessary to place more quantitative constraints on the 
spatial distribution of hot plasma and the degree of wind 
attenuation in the observed O stars; but some initial qual- 
itative inferences are already possible. For some cases, 
the implied minimal wind attenuation combined with some 
significant bulk motion are puzzling, both in the context 
of coronal X-ray emission and in the context of the strong 
radiation-driven winds that are known to exist on these 
stars from UV, optical, IR, and radio observations. One 
possible solution to the wind attenuation problem might 
be found in clumping. A "porous" wind might allow many 
of the X-rays to stream through unattenuated, yet also 
lead to significant UV absorption if the interclump wind 
medium still has sufficient column density to produce opti- 
cally thick UV lines. (The UV line transition probabilities 
are significantly larger than the bound-free X-ray cross sec- 
tions.) An alternative is that the mass- loss rates of O stars 
may have been overestimated. It does not seem likely, how- 
ever, that the winds are dense and smooth but transparent 
to X-rays, as that would require levels of ionization in the 
wind that would not allow for efficient radiation-driving. 
And, if the winds are indeed optically thin, then to dis- 
criminate between coronal and wind models will require a 
careful analysis of the shifts and asymmetries observed in 
the X-ray line profiles. 

Finally, although we developed our analysis with hot 
stars in mind, it could be extended to other X-ray sources 
in which lines are emitted within a nearly symmetric ex- 
panding medium that may also have substantial contin- 
uum opacity. This could include nova winds, planetary 
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nebulae, active galactic nuclei, and X-ray binaries. 
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